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ABSTRACT 

We investigate the ability of primordial gas clouds to retain molecular hydrogen (H2) 
during the initial phase of the reionization epoch. We find that before the Stromgren 
spheres of the individual ionizing sources overlap, the UV background below the 
ionization threshold is able to penetrate large clouds and suppress their H2 abundance. 
The consequent lack of H2 cooling could prevent the collapse and fragmentation of 
clouds with virial temperatures T v i r ^ 10 4 K (or masses ^ 1O 8 M0[(1 + -^vir) / 10] — 3 / 2 ). 
This negative feedback on structure-formation arises from the very first ionizing sources, 
and precedes the feedback due to the photoionization heating. 

Subject headings: cosmology:theory-galaxies:formation-molecular processes- 
radiative transfer 

1. Introduction 

Lately, there has been a renewed interest in the reionization of the intergalactic medium 
(IGM) due to its signature on the cosmic microwave background anisotropies (Tegmark, Silk k 
Blanchard 1994; Aghanim et al. 1996; Loeb 1996). Observations of QSO spectra imply that the 
IGM has been highly ionized already by the redshift z ~ 5 (Schneider, Schmidt k Gunn 1991). 
Evidence for the presence of an ionizing background radiation with the required magnitude is 
provided by the proximity effect (Bajtlik, Duncan k Ostriker 1988; Lu, Wolfe, k Turnshek 1991; 
Bechtold 1994 and references therein). The ionizing sources could have either been stars or quasars 
(Miralda-Escude k Ostriker 1990; Cen k Ostriker 1993; Meiksinfe Madau 1993; Shapiro, Giroux 
k Babul 1994; Fukugita k Kawasaki 1994). 

In Cold Dark Matter (CDM) cosmologies the first generation of objects forms at redshifts 
10-50. Fragmentation of these objects into stars would result in early reionization. Couchman k 
Rees (1986) have argued that the heating associated with the photoionization of the IGM would 
inevitably raise the gas temperature to T ~ 10 4 K and increase the value of the cosmological Jeans 
mass well above its initial value of ~ 1O 5 M . This effect would tend to suppress further baryonic 
structure formation until the nonlinear mass scale rises above the new Jeans mass. This negative 
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feedback serves as a buffer that monitors the fraction of baryons converted into massive stars and 
the metalicity of the IGM as a result of early reionization (Fukugita & Kawasaki 1994), and it could 
in principle explain the nearly universal metalicity value of ~ 1% observed in Lyman-a absorption 
systems (Tytler k Fan 1994; Tytler et al. 1995; Cowie et al. 1995). 

The ability of primordial gas clouds with virial temperatures T v i r ~ 10 2-4 K to cool 
radiatively and fragment after they virialize depends strongly on the existence of molecular hydrogen 
in them (Haiman, Thoul & Loeb 1996, hereafter HTL96; Tegmark et al. 1996). If the existence of 
H2 is required in order to produce ionizing sources such as stars or quasars, then it is important 
to consider the negative feedback that an ionizing background has on the H2 abundance in the 
universe. As hydrogen molecules can be destroyed by photons below the Lyman limit (Stecher & 
Williams 1967) to which the universe is transparent before being ionized, it is possible that the 
feedback due to H2 destruction would precede the feedback due to photoionization heating. 

In the early phase of reionization, after the first few sources of ionizing radiation appear, 
the Stromgren spheres of the individual sources are well separated from each other. During this 
phase, the ionizing photons from each source are absorbed at the boundary of their local Stromgren 
spheres, leaving the bulk of the Universe neutral. However, the photons below the Lyman limit 
escape the Stromgren spheres and are able to dissociate H2 molecules elsewhere. Ionizing sources 
can therefore affect the H2 abundance at a large distance (Dekel & Rees 1987). 

In a previous paper (Haiman, Rees & Loeb 1996; hereafter HRL96) we studied the effect 
of quasar-like ionizing sources (with a power law spectrum) on the H2 cooling of gas clouds after 
the universe is reionized. We have found that in very dense clouds (n ^ 10 3 cm -3 ) the ionizing 
radiation enhances the formation of molecules, and results in a net cooling of the gas. Here we 
examine the effect of the background radiation on the H2 molecules at lower densities, inside clouds 
that are collapsing when most of the universe is still neutral. In contrast with HRL96, we calculate 
the spectrum of the background radiation before the reionization epoch. The processed photon 
spectrum above the Lyman limit is strongly suppressed due to HI continuum absorption throughout 
the universe. Below the Lyman limit, the spectrum is modulated due to HI line absorption. By 
estimating the required sizes for clouds to self-shield and retain their H2 molecules, we demonstrate 
that the H2 molecules inside virialized clouds are easily destroyed before the universe gets reionized. 
The photodissociation is due to photons softer than the Lyman limit; our results are therefore 
insensitive to whether the first objects emit a thermal (0-B star) or non-thermal spectrum. The 
early destruction of H2 may provide a negative feedback that suppresses collapse and fragmentation 
prior to the feedback due to photoionization heating of neutral hydrogen. These two feedback 
scenarios could differ appreciably in the redshift at which the IGM became ionized, the predicted 
number of population III stars at the end of the reionization epoch, and the scales of the aggregates 
in which they, and hence the first nucleosynthesis products, are formed. 

The paper is organized as follows. In § 2 we calculate the characteristic spectrum of the 
ionizing radiation during the early phase of reionization, when the Stromgren spheres around 
individual sources do not overlap. In § 3 we use this spectrum to examine the question whether 
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relatively weak fluxes of photons can penetrate large clouds and destroy H2 in them. Finally, § 4 
summarizes the cosmological implications of the feedback due to the destruction of H2 molecules. 

2. Spectrum of UV Background During the Early Phase of Reionization 

We consider the first stage of the reionization epoch, when most of the universe is still 
neutral and the Stromgren spheres around individual sources do not overlap. In order to calculate 
the spectrum of the UV background outside the Stromgren spheres, we make four assumptions: 
(i) the ionizing sources turn on simultaneously at the redshift z on ; (ii) the sources are uniformly 
distributed in space; (iii) the emitted spectral intensity is a power-law with an index a; and (iv) the 
source population maintains a constant comoving density. Assumption (iii) would be appropriate if 
the ionizing sources were mainly accreting black holes; our conclusions, however, depend primarily 
on the emission of photons in the 11.18-13.6 ev range, and therefore apply with little change if the 
radiation comes from and B stars. 

The spectrum of the background radiation is suppressed above the Lyman limit through 
absorption by the column of neutral hydrogen and helium between the source and the observer. 
Photons below the Lyman limit will freely escape the Stromgren spheres around their sources. 
Since H2 is dissociated by line absorption into the Lyman and Werner lines of molecular hydrogen 
between 11.18-14.67 eV, special care is needed in calculating the spectrum in this range. 

The cumulative flux (in units of erg cm -2 s _1 ster -1 Hz -1 ) observed at the frequency 
v and redshift z \> s due to sources with a comoving emission coefficient j{y,z) (in units of 
erg cm -3 s -1 ster -1 Hz -1 ) between the redshifts z \> s and z oa can be written as 



where v z = v(l + z)/(l + £ bs)> an( i T ( z ) 1S the optical depth due to neutral hydrogen and helium 
between the observer's redshift z \> s and the redshift z, 



Above the hydrogen ionization threshold, equation (1) must be evaluated numerically. 
Below this threshold, the flux is given by the integral in equation (1), with r = if we 
ignore absorption due to lines. Assuming go = 1/2 and writing the emission coefficient as 
j{y,z) = jo{ylvjj ) -Qf , where vjj is the hydrogen ionization energy, we obtain in this case 
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In reality, the spectrum below the ionization threshold is processed due to line absorption in the 
Lyman bands of neutral hydrogen. When the IGM is still neutral, the optical depth across the 
universe in the resonant Lyman lines is ~ 10 6 , and so each photon gets absorbed at the redshift 
when its frequency equals one of the Lyman lines. Hence, photons observed to have a frequency v 
at a redshift z \> s could have only been emitted from redshifts that are between the observer and 
an effective screen located some distance away. The redshift of the screen z screen depends on t'Lyi, 
the frequency of the nearest Lyman line above is, through the resonance condition 

( 1 ~t~ ^screen) /-\ 

We can incorporate the effect of line absorption by replacing the turn-on redshift z on with the 
screen redshift z screen in equation (4). This is equivalent to multiplying the flux in equation (4) by 
the modulating factor 

l-[(l + *obs)/(l + *s. ^ 5 + " 
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So far we have ignored the radiative decay of the excited hydrogen atoms. In reality, each 
absorbing screen will re-emit line photons. Lya photons are just re-cycled until they redshift out of 
resonance, because their destruction rate by the 2s^ls two-photon decay is much slower than the 
cosmological expansion [cf. equation 6.114 in Peebles (1993)]. Therefore, there is no net absorption 
at the redshift corresponding to the Lya resonance. On the other hand, Ly/3 and higher Lyman line 
photons are re-emitted as a sum of Lya and other Balmer or lower line photons. Consequently, the 
screens corresponding to Ly/3 and higher Lyman lines produce two effects: (i) they block the view 
of all sources beyond the screen redshift, and (ii) they cascade to the ground state by emitting Lya 
and some other, lower-frequency, photons. This results in a sawtooth-modulation of the spectrum 
in the frequency range from Ly/3 up to the Lyman limit, and in a discontinuous step in the spectrum 
below the Lya frequency. 

Although it is outside the range of the Lyman and Werner bands of molecular hydrogen, the 
flux just below the Lya frequency effects the dissociation rate of H~ and H^~. To calculate this flux, 
we assume that each radiative excitation to the 3p state or higher is followed instantaneously by 
a decay to the ground state with the emission of exactly one Lya photon (plus other lower energy 
photons which are not important for our analysis). In addition, we assume that the re-emitted 
Lya photons immediately shift out of resonance due to cosmological expansion, and proceed to the 
observer without further scattering. Under these assumptions, the additional flux at the frequency 
v below Lya due to Lya emission from all screens can be written as a sum of terms, each of which 
being similar to equation (4), 
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where z a is the redshift at which the absorbing screens need to be located in order that their Lya 
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and z\ is the highest redshift at which the Lyi photons absorbed at z a could have been emitted, i.e. 
the location of a Ly(i+l)-absorbing screen, 



during absorption and re-emission. 

Figure 1 shows, as an example, the processed spectrum at z \> s = 25 assuming the sources 
turn on at z on = 35. The top panel shows that above the Lyman limit the flux is suppressed up to 
the keV energy range, and the bottom panel presents the sawtooth modulation below the Lyman 
limit, as well as the discontinuous jump at the Lya frequency. 

The assumption that the UV source distribution is homogeneous remains valid as long as 
the radius around the observer out to the nearest absorbing screen is larger than the clustering 
length of the sources. Clearly the assumption breaks down for the highest transitions in the Lyman 
series. However, our calculation is concerned with the molecular Lyman or Werner lines which are 
separated from the nearest Lyman line by redshift intervals that are typically [Az/(1 + z)] i; 10~ 2 , 
i.e. much larger than the expected clustering lengths at high redshifts. 

One may wonder whether it is possible to detect the reionization spectrum directly. The 
sawtooth template presented in figure 1 has a characteristic structure that should survive until the 
present time, provided the breakthrough epoch during which the universe becomes transparent 
to lines is sudden. Following an abrupt breakthrough epoch (as expected for an unclustered 
source population), the spectrum below the Lyman limit would be preserved and would show 
up as a slight modulation on top of the contribution from more recent sources which dominate 
the observed infrared background. The frequencies at which this sawtooth modulation appears 
depends on the reionization redshift, and its detection would provide invaluable information about 
the reionization epoch. Unfortunately, the amplitude of the sawtooth modulation is small. We 
estimate this amplitude based on the DIRBE observations (see Bond 1995 for a review), and 
assuming a constant comoving density of ionizing sources that ionize the universe abruptly at 
z = 5 and produce J 21 = 0.1 at z = 2. For these conditions the modulation amplitude is smaller 
than 10~ 6 , implying that direct detection of the reionization signal is not feasible. 
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3. Destruction of Molecular Hydrogen 



We may now use the spectrum illustrated in Fig. 1 to derive the equilibrium H2 fraction 
inside overdense regions in the universe. In particular, we would like to find out how massive should 
a gas cloud be in order to form sufficient H2 to cool on its dynamical time-scale. 
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In calculating the H2 fraction, we consider a homogeneous cloud of primordial H/He ratio 
with a density p and temperature T, which is composed of the nine species H, H~, H + , He, 
He + , He ++ , H2, H^, and e~. The chemical composition of the gas within an overdense region is 
determined by the interaction between these species and the background radiation. For the present 
purpose, we assume chemical and ionization equilibrium, which is reasonable for the densities and 
temperatures of interest. The equilibrium H2 fraction is expected to provide an upper limit on 
the actual H2 fraction and so our constraints on the destruction of H2 are conservative. The 
assumption of equilibrium is powerful in that it makes our conclusions independent of the details 
of hydrodynamic and chemical history of the gas. 

As emphasized before, the most important element of the chemical network in the present 
context is the dissociation of H2 by photons with energies below 13.6eV. The significance of this 
two-step process was first highlighted by Solomon in 1965 (cf. Field et al. 1966) and subsequently 
studied quantitatively by Stecher & Williams (1967). Here we calculate the dissociation rate for this 
process by summing the oscillator strengths for the Lyman and Werner bands of H2 from Allison 
and Dalgarno (1970), and multiplying it with the dissociation fractions (for the second step in the 
two-step process) from Dalgarno and Stephens (1970), and with the flux at each molecular line 3 . 
HTL96 have shown that the collapse of a primordial gas cloud beyond virialization depends on its 
ability to form H2 and cool radiatively. The size of a cloud that is able to form H2 depends, in 
turn, on its self- shielding in the Lyman and Werner lines. In order to correct the dissociation rate 
for self- shielding as a function of H2 column density, we used the self- shielding function from de 
Jong, Dalgarno & Boland (1980). Other details of the chemistry, including the full list of reactions 
and rate coefficients that we used, are described in HRL96. 

For the power-law spectrum in equations (4) and (7), we choose a power-law index a = 0.7, 
and an upper cutoff of E = 40keV for the background flux, as expected from quasar spectra and 
fits to the X-ray background (Fabian & Barcons 1992). We leave the normalization of the flux as 
a free parameter, and specify it as J 21, the flux (without absorption) at the Lyman limit in units 
of 10~ 21 erg s -1 cm -2 Hz -1 sr -1 . 

We calculate the equilibrium fraction of H2 within a UV-irradiated homogeneous gas cloud 
by integrating the chemical rate equations. In order to construct the equilibrium H2 profile as 
a function of radius, we subdivide the homogeneous cloud into a large number of thin plane- 
parallel slabs, and in each slab we integrate the rate equations until equilibrium is reached, taking 
account of the cumulative optical depth due to the H, He, and the H2 lines. We consider a top-hat 
overdensity of 200 relative to the background universe, and five different combinations of redshift 
and temperature. As a plausible example, we first set z=2Q and T = 10 3 K. We then repeat the 
calculation for z=10 and 50 keeping T = 10 3 K, and for T = 10 2 - 2 and 10 3 - 5 K keeping z=20. We 



3 The effect of the sawtooth modufation on the H2 dissociation rate can be fairfy farge. tn the fimit z b s /zon — ► 0, 
the modufation decreases the totaf dissociation rate by a factor of 32.5. This is the maximum effect the sawtooth 
modufation can introduce, and as z b s approaches z on , the decrease in the rate is smaffer [cf. eq. (6)]. 



- 7 - 



selected these temperatures to bracket the range favorable for H2 cooling; H2 is not formed below 
10 2 ' 2 K, while it is collisionally dissociated above ~ 10 3 ' 5 K (see HRL96). The background universe 
is assumed to have fib = 0.1 and Ho = 50 km s _1 Mpc -1 . As both the temperature and the density 
of the gas are fixed, the only free parameter left in this calculation is the normalization of the flux, 
J21. 

An example for a typical H2 fraction profile in a cloud is shown in figure 2. This figure shows 
the variation of the equilibrium H2 fraction /h2 as a function of neutral hydrogen column density Ah 
for a hydrogen density ran = 1 cm -3 , a temperature T = 1000K, and a flux J 21 = 0.1. Based on the 
processes which determine the H2 fraction, three characteristic regimes may be distinguished. Up to 
hydrogen column densities of Ah ~ 10 19 cm -2 , the most important mechanism for the destruction 
of H2 is direct photoionization by photons with energies above 15.4eV. At larger Ah, the destruction 
rate is dominated by the Solomon process, and /h2 drops. At still larger depths, Ah ~ 10 23 cm -2 , 
the Lyman and Werner lines of H2 become self- shielding, and collisional destruction by fast electrons 
(supplied mainly by the ionization of helium) dominates. The equilibrium state deep inside the 
cloud where Ah ~ 10 24 cm -2 is fully molecular, but this state is not likely to be reached within a 
Hubble time (see discussion in HRL96). 

Figure 2 also shows the ratio of the cooling time to the dynamical time i C oolAdyn as a 
function of depth, where 
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The total cooling function A tot includes photo-ionization heating, H2 cooling, and other terms that 
are summarized in detail in HRL96. The H2 cooling term was taken from Lepp & Shull 1984. 
Up to a column density Ah ~ 10 22 cm -2 , photo-ionization heating dominates H2 cooling, and the 
total cooling function, as well as t coo \ are negative, as shown in figure 2. At Ah ~ 10 22 cm -2 , 
the amount of H2 is sufficiently large so that H2 cooling is effective, and at this depth i C ool/^dyn 
becomes positive and quickly drops below unity. 

Using profiles similar to figure 2 for each value of /21 , we have found the depth rH2 at which 
a sufficient H2 abundance is formed so that the H2 cooling time is shorter than the dynamical 
time. For each value of J21, we then define the minimum cooling mass to be Mh2 = 47r/3 / or H2 . 
Figure 3a shows the cooling mass as a function of J 21 for five different combinations of redshift 
and temperature. The most probable combination T = 10 3 K, z = 20 is shown as a thick line. In 
figure 3b we show the same curves, obtained by calculating the minimum cooling mass without the 
sawtooth modulation below 13.6eV. 

Figures 3a and 3b bracket the range of possible H2 photodissociation rates, and imply that 
in general the cooling masses need to be exceedingly large, even when J 21 has a relatively low value. 
In particular, the proximity effect confirms the existence of a background flux of at least J 21 ~ 0.1 
that keeps the IGM fully ionized at z ~ 5. A similar, or even higher value of J 21 is expected at 
the time of breakthrough during the reionization epoch at higher redshifts, since the density of 
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the IGM is then higher. Figures 3a and 3b show that even at fluxes that would be able to ionize 
the IGM only partially, J 21 ^ 0.01, the mass required to shield H2 against the radiation is large, 
M m k, 1O 9 M . This conclusion can be inferred from both figures 3a and 3b, implying that it is, 
in fact, insensitive to the details of the sawtooth-modulation of the spectrum. 

4. Cosmological Implications 

Our main result is that the mass of a virialized primordial cloud has to be exceedingly high 
in order to self-shield against even a low level of background radiation. Figures 3a and b show the 
mass required to form enough H2 molecules to cool the cloud on its dynamical time-scale, versus 
the normalization of the flux J 21 ■ The minimum cooling masses shown in these figures are indeed 
high, in all cases they are above ~ 10 9 M already for fluxes that are too low to fully ionize the 
IGM. 

It has been argued in previous works (HTL96; Tegmark et al. 1996) that clouds with virial 
temperatures 100 K ^ T v - ir ^ 10 4 K are able to cool by H2 and continue to collapse after they 
virialize. These papers, however, did not take into account any external UV radiation. Our present 
results imply that the appearance of the very first UV sources will suppress the H2 abundance in 
these objects, and prevent their post-virialization collapse. 

The cooling masses we obtained depend primarily on the number of photons in the narrow 
spectral range of 11.18eV-13.6eV. Our conclusions are therefore insensitive to the detailed shape 
of the spectrum, and should hold for either stellar (blackbody) or quasar (power-law) spectra. We 
also note that in our calculation of the cooling mass, we assumed a homogeneous spherical cloud. 
Although centrally condensed clouds are, in general, more effective than homogeneous clouds in 
self- shielding their molecules, a steepening of the density profile is unlikely to decrease the value of 
the required cooling mass by orders of magnitude. 

We have presented our results for a range of temperatures and redshifts. Our results are 
relevant to any model where the initial fluctuations have amplitudes decreasing with scale, so that 
cosmic structures form hierarchically. Such models predict different turn-on redshifts for the first 
sources of radiation; in PIB models, turn-on may occur at z > 100; in CDM in the range z = 10 — 50; 
for mixed dark matter models the first structures may form still more recently. Molecular cooling 
is more efficient at high densities, and therefore at larger redshifts; however in all these models 
it determines the scale of the first objects that condense out and contribute the first injection of 
heat into the universe. The amount of UV generated by each such object is very uncertain; it 
depends on the efficiency of star formation, and whether the IMF favors massive stars, or even 
super-massive objects or black holes. However, we have shown that H2 cooling is quenched even in 
massive objects before the UV background has reached the level needed to photoionize the entire 
IGM. Therefore the radiation of the first sources exerts an important feedback on the cosmogonic 
process: the only clouds able to fragment into new stars are those with virial temperature above 
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~ 10 4 K, i.e. with masses ^ 1O 8 M0[(1 + -^vir) / 10] — 3 / 2 ), so that isothermal contraction is possible 
due to atomic line cooling. Clouds that virialize at lower temperatures will, in the absence of a 
sufficient amount of H2, maintain virial equilibrium and not fragment into stars. 

We therefore draw the robust conclusion that the IGM remained predominantly neutral 
until a sufficient number of objects with virial temperatures above fO 4 K (masses above ~ 
fO 8 M [(f + ^ vir )/f0]" 3 / 2 ) had gone non-linear. Most of the 0-B stars or accreting black holes 
that photoionized the IGM had to form in systems at least as large as this. Formation of such 
systems would have continued unimpeded until ionization was complete; subsequently J 21 rises 
sharply to a value of order unity, when the universe becomes, in effect, a single HII region. The 
only net cooling of a fully photoionized gas comes from bremsstrahlung, which is less effective than 
the collisionally-excited line emission from a partially ionized gas (Efstathiou, 1992; Weinberg, 
Hernquist and Katz 1996, Navarro and Steinmetz 1996). The completion of photoionization may 
therefore signal another pause in the cosmogonic process, associated with the further decline in the 
efficiency of cooling, and the consequent further increase in the minimum scale that can collapse. 

We thank Alex Dalgarno for discussions on the chemistry, and the referee, Alain Blanchard, 
for useful comments. AL acknowledges support from the NASA ATP grant NAG5-3085. 
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Fig. 1. — The average flux (in units of erg cm -2 s _1 ster -1 Hz -1 ) during the initial phase of 
the reionization epoch. The upper panel shows that absorption by neutral hydrogen and helium 
suppresses the flux above 13.6eV up to the keV range. The lower panel shows the sawtooth- 
modulation due to line absorption below 13.6eV. 
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Fig. 2. — The equilibrium profile of the H2 fraction fu 2 , and the ratio of the cooling time to the 
dynamical time i C ool/^dyn- At low column densities, photoionization results in a net heating, this 
is expressed by t coo \ < 0. At the column density Ah ~ 10 21 cm -2 the H2 Lyman and Werner 
lines become self- shielding and the molecular fraction starts to rise. When the molecular fraction 
increases above 10~ 4 , H2 cooling dominates over photoionization heating, resulting in net cooling 
(^cool > 0). At a slightly larger depth, the H2 cooling time decreases below the dynamical time. 
The column density at which this occurs, Ah ~ 10 22 cm -2 , defines the minimum cooling mass used 
in § 3. 
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Fig. 3. — (a) The minimum mass Mh2 of a homogeneous, spherically symmetric cloud which is 
able to self-shield and form sufficient amount of H2 at its center so that it cools on a dynamical 
time-scale. The lines show Mh2 for five different combinations of redshift and temperature. Among 
these, the intermediate combination (log T,z) = (3,20) is shown as a thick line, (b) Same as (a), 
but with the minimum mass calculated without the sawtooth-modulation of the spectrum below 
13.6eV. Note that the cooling masses in (a) and (b) differ only slightly. In both cases, they are 
high already at fluxes that are too low to ionize the IGM. 



